Detecting chromosomal interactions in Capture Hi-C data with CHiCAGO and companion tools.

Capture Hi-C is widely used to obtain high-resolution profiles of chromosomal interactions involving, at least on one end, regions of interest such as gene promoters. Signal detection in Capture Hi-C data is challenging and cannot be adequately accomplished with tools developed for other chromosome conformation capture methods, including standard Hi-C. Capture Hi-C Analysis of Genomic Organization (CHiCAGO) is a computational pipeline developed specifically for Capture Hi-C analysis. It implements a statistical model accounting for biological and technical background components, as well as bespoke normalization and multiple testing procedures for this data type. Here we provide a step-by-step guide to the CHiCAGO workflow that is aimed at users with basic experience of the command line and R. We also describe more advanced strategies for tuning the key parameters for custom experiments and provide guidance on data preprocessing and downstream analysis using companion tools. In a typical experiment, CHiCAGO takes ~2-3 h to run, although pre- and postprocessing steps may take much longer.

An early cell shape transition drives evolutionary expansion of the human forebrain.

The human brain has undergone rapid expansion since humans diverged from other great apes, but the mechanism of this human-specific enlargement is still unknown. Here, we use cerebral organoids derived from human, gorilla, and chimpanzee cells to study developmental mechanisms driving evolutionary brain expansion. We find that neuroepithelial differentiation is a protracted process in apes, involving a previously unrecognized transition state characterized by a change in cell shape. Furthermore, we show that human organoids are larger due to a delay in this transition, associated with differences in interkinetic nuclear migration and cell cycle length. Comparative RNA sequencing (RNA-seq) reveals differences in expression dynamics of cell morphogenesis factors, including ZEB2, a known epithelial-mesenchymal transition regulator. We show that ZEB2 promotes neuroepithelial transition, and its manipulation and downstream signaling leads to acquisition of nonhuman ape architecture in the human context and vice versa, establishing an important role for neuroepithelial cell shape in human brain expansion.

Ribosome collisions trigger cis-acting feedback inhibition of translation initiation.

Translation of aberrant mRNAs can cause ribosomes to stall, leading to collisions with trailing ribosomes. Collided ribosomes are specifically recognised by ZNF598 to initiate protein and mRNA quality control pathways. Here we found using quantitative proteomics of collided ribosomes that EDF1 is a ZNF598-independent sensor of ribosome collisions. EDF1 stabilises GIGYF2 at collisions to inhibit translation initiation in cis via 4EHP. The GIGYF2 axis acts independently of the ZNF598 axis, but each pathway's output is more pronounced without the other. We propose that the widely conserved and highly abundant EDF1 monitors the transcriptome for excessive ribosome density, then triggers a GIGYF2-mediated response to locally and temporarily reduce ribosome loading. Only when collisions persist is translation abandoned to initiate ZNF598-dependent quality control. This tiered response to ribosome collisions would allow cells to dynamically tune translation rates while ensuring fidelity of the resulting protein products.

Genetically Defined Functional Modules for Spatial Orienting in the Mouse Superior Colliculus.

In order to explore and interact with their surroundings, animals need to orient toward specific positions in space. Throughout the animal kingdom, head movements represent a primary form of orienting behavior. The superior colliculus (SC) is a fundamental structure for the generation of orienting responses, but how genetically distinct groups of collicular neurons contribute to these spatially tuned behaviors remains largely to be defined. Here, through the genetic dissection of the murine SC, we identify a functionally and genetically homogeneous subclass of glutamatergic neurons defined by the expression of the paired-like homeodomain transcription factor Pitx2. We show that the optogenetic stimulation of Pitx2ON neurons drives three-dimensional head displacements characterized by stepwise, saccade-like kinematics. Furthermore, during naturalistic foraging behavior, the activity of Pitx2ON neurons precedes and predicts the onset of spatially tuned head movements. Intriguingly, we reveal that Pitx2ON neurons are clustered in an orderly array of anatomical modules that tile the entire intermediate layer of the SC. Such a modular organization gives origin to a discrete and discontinuous representation of the motor space, with each Pitx2ON module subtending a defined portion of the animal's egocentric space. The modularity of Pitx2ON neurons provides an anatomical substrate for the convergence of spatially coherent sensory and motor signals of cortical and subcortical origins, thereby promoting the recruitment of appropriate movement vectors. Overall, these data support the view of the superior colliculus as a selectively addressable and modularly organized spatial-motor register.

Epigenetic Regulation of Vascular Smooth Muscle Cells by Histone H3 Lysine 9 Dimethylation Attenuates Target Gene-Induction by Inflammatory Signaling.

OBJECTIVE: Vascular inflammation underlies cardiovascular disease. Vascular smooth muscle cells (VSMCs) upregulate selective genes, including MMPs (matrix metalloproteinases) and proinflammatory cytokines upon local inflammation, which directly contribute to vascular disease and adverse clinical outcome. Identification of factors controlling VSMC responses to inflammation is therefore of considerable therapeutic importance. Here, we determine the role of Histone H3 lysine 9 di-methylation (H3K9me2), a repressive epigenetic mark that is reduced in atherosclerotic lesions, in regulating the VSMC inflammatory response. Approach and Results: We used VSMC-lineage tracing to reveal reduced H3K9me2 levels in VSMCs of arteries after injury and in atherosclerotic lesions compared with control vessels. Intriguingly, chromatin immunoprecipitation showed H3K9me2 enrichment at a subset of inflammation-responsive gene promoters, including MMP3, MMP9, MMP12, and IL6, in mouse and human VSMCs. Inhibition of G9A/GLP (G9A-like protein), the primary enzymes responsible for H3K9me2, significantly potentiated inflammation-induced gene induction in vitro and in vivo without altering NFκB (nuclear factor kappa-light-chain-enhancer of activated B cell) and MAPK (mitogen-activated protein kinase) signaling. Rather, reduced G9A/GLP activity enhanced inflammation-induced binding of transcription factors NFκB-p65 and cJUN to H3K9me2 target gene promoters MMP3 and IL6. Taken together, these results suggest that promoter-associated H3K9me2 directly attenuates the induction of target genes in response to inflammation in human VSMCs. CONCLUSIONS: This study implicates H3K9me2 in regulating the proinflammatory VSMC phenotype. Our findings suggest that reduced H3K9me2 in disease enhance binding of NFκB and AP-1 (activator protein-1) transcription factors at specific inflammation-responsive genes to augment proinflammatory stimuli in VSMC. Therefore, H3K9me2-regulation could be targeted clinically to limit expression of MMPs and IL6, which are induced in vascular disease.

Author Correction: Promoter interactome of human embryonic stem cell-derived cardiomyocytes connects GWAS regions to cardiac gene networks.

In the original version of the Article, the gene symbol for tissue factor pathway inhibitor was inadvertently given as 'TFP1' instead of 'TFPI'. This has now been corrected in both the PDF and HTML versions of the Article.

Genome organization and chromatin analysis identify transcriptional downregulation of insulin-like growth factor signaling as a hallmark of aging in developing B cells.

BACKGROUND: Aging is characterized by loss of function of the adaptive immune system, but the underlying causes are poorly understood. To assess the molecular effects of aging on B cell development, we profiled gene expression and chromatin features genome-wide, including histone modifications and chromosome conformation, in bone marrow pro-B and pre-B cells from young and aged mice. RESULTS: Our analysis reveals that the expression levels of most genes are generally preserved in B cell precursors isolated from aged compared with young mice. Nonetheless, age-specific expression changes are observed at numerous genes, including microRNA encoding genes. Importantly, these changes are underpinned by multi-layered alterations in chromatin structure, including chromatin accessibility, histone modifications, long-range promoter interactions, and nuclear compartmentalization. Previous work has shown that differentiation is linked to changes in promoter-regulatory element interactions. We find that aging in B cell precursors is accompanied by rewiring of such interactions. We identify transcriptional downregulation of components of the insulin-like growth factor signaling pathway, in particular downregulation of Irs1 and upregulation of Let-7 microRNA expression, as a signature of the aged phenotype. These changes in expression are associated with specific alterations in H3K27me3 occupancy, suggesting that Polycomb-mediated repression plays a role in precursor B cell aging. CONCLUSIONS: Changes in chromatin and 3D genome organization play an important role in shaping the altered gene expression profile of aged precursor B cells. Components of the insulin-like growth factor signaling pathways are key targets of epigenetic regulation in aging in bone marrow B cell precursors.

Promoter interactome of human embryonic stem cell-derived cardiomyocytes connects GWAS regions to cardiac gene networks.

Long-range chromosomal interactions bring distal regulatory elements and promoters together to regulate gene expression in biological processes. By performing promoter capture Hi-C (PCHi-C) on human embryonic stem cell-derived cardiomyocytes (hESC-CMs), we show that such promoter interactions are a key mechanism by which enhancers contact their target genes after hESC-CM differentiation from hESCs. We also show that the promoter interactome of hESC-CMs is associated with expression quantitative trait loci (eQTLs) in cardiac left ventricular tissue; captures the dynamic process of genome reorganisation after hESC-CM differentiation; overlaps genome-wide association study (GWAS) regions associated with heart rate; and identifies new candidate genes in such regions. These findings indicate that regulatory elements in hESC-CMs identified by our approach control gene expression involved in ventricular conduction and rhythm of the heart. The study of promoter interactions in other hESC-derived cell types may be of utility in functional investigation of GWAS-associated regions.

Long-Range Enhancer Interactions Are Prevalent in Mouse Embryonic Stem Cells and Are Reorganized upon Pluripotent State Transition.

Transcriptional enhancers, including super-enhancers (SEs), form physical interactions with promoters to regulate cell-type-specific gene expression. SEs are characterized by high transcription factor occupancy and large domains of active chromatin, and they are commonly assigned to target promoters using computational predictions. How promoter-SE interactions change upon cell state transitions, and whether transcription factors maintain SE interactions, have not been reported. Here, we used promoter-capture Hi-C to identify promoters that interact with SEs in mouse embryonic stem cells (ESCs). We found that SEs form complex, spatial networks in which individual SEs contact multiple promoters, and a rewiring of promoter-SE interactions occurs between pluripotent states. We also show that long-range promoter-SE interactions are more prevalent in ESCs than in epiblast stem cells (EpiSCs) or Nanog-deficient ESCs. We conclude that SEs form cell-type-specific interaction networks that are partly dependent on core transcription factors, thereby providing insights into the gene regulatory organization of pluripotent cells.

Global reorganisation of cis-regulatory units upon lineage commitment of human embryonic stem cells.

Long-range cis-regulatory elements such as enhancers coordinate cell-specific transcriptional programmes by engaging in DNA looping interactions with target promoters. Deciphering the interplay between the promoter connectivity and activity of cis-regulatory elements during lineage commitment is crucial for understanding developmental transcriptional control. Here, we use Promoter Capture Hi-C to generate a high-resolution atlas of chromosomal interactions involving ~22,000 gene promoters in human pluripotent and lineage-committed cells, identifying putative target genes for known and predicted enhancer elements. We reveal extensive dynamics of cis-regulatory contacts upon lineage commitment, including the acquisition and loss of promoter interactions. This spatial rewiring occurs preferentially with predicted changes in the activity of cis-regulatory elements and is associated with changes in target gene expression. Our results provide a global and integrated view of promoter interactome dynamics during lineage commitment of human pluripotent cells.

Lineage-Specific Genome Architecture Links Enhancers and Non-coding Disease Variants to Target Gene Promoters.

Long-range interactions between regulatory elements and gene promoters play key roles in transcriptional regulation. The vast majority of interactions are uncharted, constituting a major missing link in understanding genome control. Here, we use promoter capture Hi-C to identify interacting regions of 31,253 promoters in 17 human primary hematopoietic cell types. We show that promoter interactions are highly cell type specific and enriched for links between active promoters and epigenetically marked enhancers. Promoter interactomes reflect lineage relationships of the hematopoietic tree, consistent with dynamic remodeling of nuclear architecture during differentiation. Interacting regions are enriched in genetic variants linked with altered expression of genes they contact, highlighting their functional role. We exploit this rich resource to connect non-coding disease variants to putative target promoters, prioritizing thousands of disease-candidate genes and implicating disease pathways. Our results demonstrate the power of primary cell promoter interactomes to reveal insights into genomic regulatory mechanisms underlying common diseases.

CHiCAGO: robust detection of DNA looping interactions in Capture Hi-C data.

Capture Hi-C (CHi-C) is a method for profiling chromosomal interactions involving targeted regions of interest, such as gene promoters, globally and at high resolution. Signal detection in CHi-C data involves a number of statistical challenges that are not observed when using other Hi-C-like techniques. We present a background model and algorithms for normalisation and multiple testing that are specifically adapted to CHi-C experiments. We implement these procedures in CHiCAGO ( http://regulatorygenomicsgroup.org/chicago ), an open-source package for robust interaction detection in CHi-C. We validate CHiCAGO by showing that promoter-interacting regions detected with this method are enriched for regulatory features and disease-associated SNPs.